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o : Abstract 

We present a global study of the simplest scalar phantom dark matter model. The best fit 
parameters of the model are determined by simultaneously imposing (i) relic density constraint 
from WMAP, (ii) 225 live days data from direct experiment XENON100, (iii) upper limit of 
gamma-ray flux from Fermi-LAT indirect detection based on dwarf spheroidal satellite galaxies, 
and (iv) the Higgs boson candidate with a mass about 125 GeV and its invisible branching ratio no 
larger than 40% if the decay of the Higgs boson into a pair of dark matter is kinematically allowed. 
The allowed parameter space is then used to predict annihilation cross sections for gamma-ray 
lines, event rates for three processes mono-6 jet, single charged lepton and two charged leptons 
plus missing energies at the Large Hadron Collider, as well as to evaluate the muon anomalous 
magnetic dipole moment for the model. 
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I. INTRODUCTION 



Evidences for the existence of dark matter are mainly coming from cosmological obser- 
vations related to the physics of gravity. These include the relic density of dark matter, 
anisotropies in the Cosmic Microwave Background (CMB), large scale structure of the uni- 
verse, as well as the bullet clusters and the associated gravitational lensing effects. While 
we still do not know what the nature of dark matter is, it is clear that there is no room to 
accommodate dark matter in the standard model (SM) of particle physics based on gauge 
invariance of SU(3)c x SU (2)^ x U (l)y and Einstein-Hilbert gravity theory based on general 
coordinate invariance. While it is plausible that the nature of dark matter may have a purely 
gravitational origin, theories that have been put forward thus far are not as convincing as 
those from the particle physics point of view. In particular the relic density strongly suggests 
that dark matter may be a weakly interacting massive particle (WIMP). If dark matter can 
indeed be related to weak scale physics, there may be hope for us to detect them in various 
underground experiments of direct detection as well as in space experiments using balloons, 
satellites, or space station of indirect detection. Furthermore, WIMP dark matter might 
be produced directly at the Large Hadron Collider (LHC) by manifesting itself as missing 
energy with a spectrum that may be discriminated from standard model background of 
neutrinos. 

In this paper, we will focus on the simplest dark matter model []] which is based on 
adding a real singlet scalar field to the SM. The communication between the scalar dark 
matter (dubbed scalar phantom in and the SM gauge bosons and fermions must then go 
through the SM Higgs boson. While there have been many studies for this simple model and 
its variants in the literature j^-Q], we believe a global study of this model is still missing. 
In this work, we will fill this gap. We use the current experimental constraints of relic 



density from WMAP js], 2 



9|, diffuse 



gamma-ray flux from indirect detection experiment of Fermi-LAT using the dwarf spheroidal 



satellite galaxies (dSphs) 10|, 111 L and a Higgs boson candidate with mass about 125 GeV 



reported recently by the LHC 12J, ll3j to deduce the best fit parameters of the model. The 
deduced parameters are used to predict various phenomenology of the model at the LHC, 
including production of the mono-6 jet, single charged lepton, and two charged leptons plus 
missing energies. We also evaluate the muon anomalous magnetic dipole moment which is 



a two loop process in the model. For a global fitting based on effective operators approach, 
see our recent work in fl4l ]. A similar global analysis for isospin violating dark matter is 
presented in [15|. 

In the next section, we will briefly review the scalar phantom model of dark matter. In 
section III, we present the global fitting for the relevant parameters of the model using the 
various experimental constraints described above. In section IV, we discuss collider phe- 
nomenology and the muon anomalous magnetic dipole moment of the model. We conclude 
in section V. Some analytical formulas of the matrix elements needed in our analysis as 
well as the expression for the muon anomalous magnetic dipole moment are collected in the 
Appendix. 



II. THE SCALAR PHANTOM MODEL 

The simplest dark matter model (SZ) [l| (dubbed scalar phantom by the authors in 
is obtained by adding one real singlet scalar x i n addition to the Higgs doublet $ to the SM. 
The scalar part of the Lagrangian is given by 

£ scalar = (D^ (D,<S>) - A (V$ - + IpxdpX ~ \™\ 2 ~ Ivx 4 ~ \px 2 &® • (1) 

A discrete Z2 symmetry of \ — > —\ while keeping all SM fields unchanged has been imposed 
to eliminate the x-, X^^> an d X 3 terms. As a result it guarantees the stability of the x 
particle and hence it may be a viable candidate for WIMP (weakly interacting massive 
particle) dark matter. Note that the x 4 term in Eq.flT]) implies a contact interaction vertex 
among the scalar dark matter. 

The virtue of this model is its simplicity. Indeed, it represents the simplest realization of a 
broad class of models, in which we could add any number of singlet scalar x to the standard 
model, or the standard model augmented by a private Higgs sector The analysis given 
here is in the spirit of seeing whether or not the simplest version of this kind of model could 
now be ruled out. 

After electroweak symmetry breaking, $ develops a vacuum expectation value v/\/2, 
where v = fi/VX = 246 GeV. After making the shift $(x) T = (0 , v + H{xj) /a/2, the 
physical Higgs field H obtains a mass m H = \/2Xv = y/2/j, and the last term in Eq.([T]) 
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becomes 



1 2 2 1 ej2 1 rr2 2 



(2) 



The first term on the right handed side of Eq.(j2]) implies the dark matter \ a ls° pick U P an 
additional contribution of \pv 2 to its mass, thus mi = m 2 + \pv 2 . We will assume m 2 x is 
always positive so that the Z 2 symmetry will never be broken, except perhaps due to black 
hole effects. The second term in Eq. (T2J) tells us that the dark matter \ can communicate 
to the SM fields and self-interact with itself via a tree level Higgs exchange, while the last 
term contributes to the relic density calculation from the process XX ~~ HH if kinematically 
allowed. If kinematics permits, the second term also allows Higgs boson to decay into a pair 
of giving rise to the invisible Higgs width. Implication of invisible Higgs width in the 
Higgs search at the LHC will be discussed further in the following sections. 

There are a few theoretical restrictions on the model, including vacuum stability, unitarity, 
and triviality. Stability of the vacuum requires the scalar potential be bounded from below. 
At tree level, we have 



2 2A?7 

A>0,r/>0,p 2 < — - 



(3) 



Tree level perturbative unitarity constraints can be deduced by considering the longitudinal 
gauge boson scatterings 17] as well as all scalar-scalar scatterings j^J 



ml < ?V ~ (712 GeVf 



T] < 8ir , and \p\ < 8n . 



(4) 



BQ. 



Analysis of the triviality of this model can be found in the literature 

to resolve some conflicts between 



Self-interacting cold dark matter was proposed in 
actual observations and WIMP theory which predicts overly dense cores in the center of 
galaxies and clusters and an overly large number of halos within the Local Group. The 



Spergel-Steinhardt bound 18| for collisional (self-interacting) dark matter, 



2 x 10 3 GeV" 3 < < 3 x 1Q 4 GeV - 3 



m. 



(5) 



can be used to constrain the contact self-coupling 77 of the scalar_p 
coupling p through Higgs exchange). We refer to previous works 



rantom (as well as the 
191 ] on this issue. 



III. GLOBAL FITTING 



In this section, we consider the global constraints coming from WMAP relic density |8|, 



the XENON100 data 



9], the Fermi- L AT upper limit o: 



diffuse gamma-ray flux based on 



dSphs (lOj, the LHC Higgs mass around 125 GeV [12|, [13|, and the upper limit of Higgs 
invisible decay branching ratio 20|. The relevant parameters of the model are the dark 
matter mass m x , the coupling p and the Higgs mass run, which the likelihood functions will 
depend on. The publicly available software package MicrOMEGAs version 2.4.5 2l| is used 
to calculate the relic density and gamma-ray flux. 



A. WMAP Relic Density 

In the calculation of the relic density, we consider all of the following two body tree level 
processes 

XX # W ^ //, W + W~, ZZ, HE . (6) 



The relic density from the WMAP 7-year result js| is 

Q c h 2 = 0.112 ±0.0056 , (7) 

where Q c is the density of the cold dark matter normalized to the critical density and h is the 
Hubble rate in unit of 100 km sec -1 Mpc" 1 . In the fitting, we use the Gaussian distribution 
for the WMAP relic density in our likelihood function, 

£relic = e 2 , (8) 

with the x 2 defined as Q 

2 (prediction — experimental central value) 2 

x = > ( 9 ) 

where the a can be read off from the WMAP experimental error jg|. We also assume a 
theoretical uncertainty r to be 10% of the prediction in order to account for the discrepancy 
due to different methods being used to solve the Boltzmann equation in different relic density 
computation packages. 



1 We apologize for abusing our notation of \ to stand for the scalar phantom dark matter as well as to 
define the chi-squared. 
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B. Direct Detection 



Here we need to calculate the elastic cross section for 

X M — > xA/- , (10) 

where M is a nucleus with Z protons and {A — Z) neutrons. For a Higgs mass around 125 
GeV, it is heavy enough to be integrated out to give an effective local interaction between the 
dark matter and the quarks. Since the local velocity of the dark matter is about v x ~ lCT 3 c, 
non-relativistic reduction is appropriate. For the present model there is no spin-dependent 
cross section for the above elastic scattering, because the dark matter is a scalar particle. 
The spin-independent cross section at zero recoil energy can be obtained as 



°Sr(0) = ^U z fp + ( A ~ z )fn\ 2 , (11) 



where fi x jy = m x mj^/( y m x + mj\f) is the reduced mass for the xN system and 

h-P^ri £ C + ^M (12) 

*■ n {q=u,d,s J 

with uq denotes the number of heavy quarks and N = p or n. The factors and f^g are 



hadronic matrix elements, and we will use their default values given in MicrOMEGAs [21]. 

For direct detection we choose the truncated Maxwell velocity distribution, which is the 
default choice in MicrOMEGAs. Because XENON100 is a counting experiment, the best 
choice for the likelihood function is Poisson distribution, 

r ,r e ~ (s+b) ( S + V° run 

•'-'direct OC • 

The b = 1.0 and o = 2.0 are the number of background and observation events taken from 



XENON100 3, respectively. In order to achieve the minimum of Xdirect equals to zero, we 
normalized Eq. f|T3|) by the factor e~° (o)° jo\. For simplicity we do not take the background 
uncertainties into account in our analysis. The signal s equals e x A r (unbiased) in the nuclear 
recoil energy range of 6.6 - 30.5 keV nr , where e is the detector efficiency once the experimental 
cuts are applied to the total number of unbiased events A^(unbiased). In other words, e is 
the fraction of A r (unbiased) generated by the Monte-Carlo simulation which survives the 
various cuts taken in XENON 100. 
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We use the following approximate method to evaluate the efficiency e as a function of 
m x . First, from XENON100 data of b = 1.0 and o = 2.0 at 90%C.L. one can obtain a 



reference s* = 5.72 using the likelihood function just mentioned. Then, taking the values 



of m x and a^ p along the 90% CL. line from Fig. 3 in (9] as input to MicrOMEGAs, we 
can obtain iV(unbiased) which is a function of m x and a^L. Hence, assuming the efficiency 
e depends only on m x and a^, it is simply given by the ratio s*/iV(unbiased). Below when 
we compute the number of events for our model, we will multiply it by this e to obtain the 
corresponding number of signals. 

For Poisson distribution we evaluate the effective \ 2 according to the following expression 

X 2 = -21n£. (14) 

C. Indirect Detection 



Since the couplings of the scalar dark matter to the SM fermions must go through the 
Higgs exchange in the SZ model, the positron flux and the anti-proton flux are probably 
too small for indirect detection! 2 . We will focus on gamma-ray flux in indirect detection \ 
For continuum gamma-ray, the relevant processes in the model are //, W + W~ , ZZ, and 
HH channels. The W + W~ , ZZ, and HH will subsequently decay into quarks and charged 
leptons which can produce gamma-rays. However, the HH final state will be discarded in 
our later numerical analysis since micrOMEGAs does not provide the Pythia decay table 



of XX ~~ ^ HH gamma-rays. Furthermore, Fermi-LAT (lOj has neither the limit 

for HH nor ZZ final state. Nevertheless, the HH contribution to gamma rays is expected 
to be small, since the Higgs decays into SM fermion pairs are suppressed by the fermion 
mass. When m x > raw, the dominant contribution comes from W + W~ , and the next 
largest contribution comes from the heaviest fermion pair that the XX can annihilate into. 
For lower DM mass, the dominant process is the following annihilation 

xx -> # w -> qq -> -> 2 7 X (15) 

where the neutral pion coming from quark fragmentation decays into two photons. 



2 However, it was shown in 



22j that antimatter signals for the present model might be promising at the 



AMS experiment on the space station. 
3 An earlier analysis of the gamma-ray signals for the present model can be found in [23| 
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Detection of one or more spectral lines would be the smoking gun signal for dark matter. 
The following annihilation processes XX ~^ 77 > 7^ > lH can give rise to discrete gamma 
lines. In the SZ model, only the following two processes 

XX -> # W -> 77 , 7^ (16) 

are possible Q. The 7Z final state is possible if m x > mz/2. The gamma-ray lines are 
located approximately at energies _E 7 ~ m x and m x (l — m 2 z /Am x ) for the two final states 
77 and 7Z, respectively, with corrections of order (v x /c) 2 ~ 10~ 6 , which is minuscule. Since 
these processes are one-loop induced and therefore suppressed, we do not include them in 
the global fitting. Instead, we will compute these cross sections after the scan and compare 
with the Fermi-LAT limits 24]. 



The dSphs of our Milky Way are satellite systems without active star formation or de- 
tected gas content. They are thus fainter and expected to be dominated by DM due to their 
own gravitational binding. Although the expected flux of gamma-rays is not as high as the 
Galactic Center, these dwarf galaxies may have a better signal-to-noise ratio. Currently, the 
most stringent upper limits on the DM annihilation cross sections in various channels are 
derived by Fermi-LAT Collaboration using the new 24 months data set with the following 



two improvements on their analysis 



10j. First, they performed a joint likelihood analysis 



to 10 satellite galaxies which can improve their statistical power. Second, they included the 
uncertainties in the dark matter distribution in these satellites entered in the astrophysical 
J factor 



j(V>) = / didn P 2 [i(^)] , (17) 

J lino— of— sight, AC 

which is the line-of-sight integral of the squared DM density, p, toward an observational 
direction, ift, integrated over a sustained solid angle, AQ. The gamma-rays flux is then 
given by 

<j>(E^) = J—^N^E)^) , (18) 

where (<jv x ) is the velocity-averaged pair annihilation cross section and N^(E) is the gamma- 
ray energy distribution per annihilation. Based on these two improvements, robust upper 
limits of 95% C.L. on the av x for the bb, t + t~ , and W + W~ channels are derived 



4 By charge conjugation, it is impossible to construct gauge invariance operators using one single photon 
field strength with arbitrary numbers of Higgs fields and partial derivatives. 
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Ill 



10j. We will use these constraints on the diffuse gamma-ray flux in our global fitting. 



Since each limit was obtained by assuming the dominance by one single channel, we can 
approximately reconstruct the upper limit suitable to our case by applying the same method 



as in Sec. 4.1 of 



25). 



In our analysis for dSphs, we adopt our likelihood function as follows 

—(^^V-fiF)- (i9) 

where erfc = 1 — erf is the complementary error function and the effective y 2 is the same 
as Eq. fll4j) . In addition, because the astrophysical J factor is expected to have a 3% un- 
certainty and the hadronization/decay tables in either MicrOMEGAs 21] or DarkSUSY 
26| have a factor of 2 uncertainty, we can then include the theoretical uncertainties as 



r = V0.03 2 + 2 2 x <7f X 95, where av x g 5 is the 95% C.L. of the reconstructed upper limit for 
our DM pair annihilation cross section. The effective x 2 f° r this case is taken to be same as 
Eq.(P3D- 



D. Higgs Mass and Its Invisible Width 

brce our scan to go to the m# ~ 125 GeV region, as suggested by the recent 



In order to 
LHC data 121 



13|, we use the Gaussian likelihood function for the Higgs mass with a central 
value of 125.3 GeV and an experimental uncertainty o ~ 0.6 GeV. Since the m# is an input, 
we do not introduce any theoretical error for the Higgs mass. Therefore, the likelihood 



function for the Higgs mass is 

£ ffig g S = e~5^^ '- . (20) 



(m H -125.3 GeV) 2 



It has been pointed out recently in [20J that the monojet search at the LHC has strongly 
disfavored B im/ = T inv / (T SM + r inv ) > 0.4 for Higgs-portal dark matter model where r SM is 
the total SM Higgs width. The invisible width of SM Higgs in the SZ model is 

1 P V / 4m 2 V 

^(H^XX) = ^-—(1--^) ■ (21) 
327T mH \ / 

We note that the invisible decay mode is not dominant in most of dark matter mass range. 
Only when m# < 130 GeV and m# > 2m x , the invisible decay of Higgs becomes significant. 
Hence, we implement the Higgs invisible decay as a 0/1 hard cut. If m# < 130 GeV, 
m H > 2m x , and B inv (H — > \X) > 0-4, we multiply £ffigg S by 0, otherwise by 1. 
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Cheung, Tseng, Tsai, Yuan, and Zee (2012) 
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FIG. 1. The profile likelihood of {m 



density [8|, XENON100 
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re glob al fitting using WMAP relic 



10[] and a 125 GeV Higgs 



12 



131 ] with an invisible branching 



E. Parameter Scan 

Engaging with MultiNest v2.7 {27I with 10000 living points, a stop tolerance factor 0.001, 
and an enlargement factor reduction parameter 0.5, we perform a random scan in the three 
dimensional parameter space of m x , p, and m# restricted in the following ranges 



1.0 < log 10 [m x / GeV] < 3.0 
-3.0 < log 10 [p] < 0.0 
114.0 GeV < m H < 130.0 GeV 



(22) 



The selected scan range of p is much smaller than the theoretical limit |p| < 8tt because 
the WMAP window is very small which only allows p < 1. Furthermore, in order to scan 
efficiently in the Higgs resonance region and cover the low p region, we use the log priors for 
m x and p as specified in Eq.( l22]) . Similar results are found for the case of negative p and 
will not be shown here. 
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After hitting the stop criteria, we collect total 440682 samples, and plot 68% and 95% 



profile likelihood confidence 
Nested Sampling algorithm 



imit contours based on 138017 samples which are selected by 
28 1 . The 68% and 95% confidence limit means that the total 



likelihood is greater than 0.32 * £(BestFit) and 0.05 * £ (Best Fit), respectively. The total 
likelihood function for our global fitting will be taken as 

£tot £relic X £direct X £indirect X £jliggs > (^^) 

and the effective total Xtot * s given by 

Xt 2 ot = -21n£ tot . (24) 

The result of the profile likelihood projected on the (m x , p) plane is shown in Fig. ([I]). 
We can clearly see that there are two branches: the vertical branch at low m x region and 
the horizontal branch hooked around at m x >100 GeV. The shape of these two branches 
is mainly due to the relic density constraint. However, XENON100 and dSphs also play a 
significant role at the junction of the two branches, p ~ 0.04 — 0.1 and 50 < m x /GeV < 200, 
where relatively large a?L and av can be easily produced. Furthermore, the hard cut due 
to the Higgs invisible branching ratio can remove some of the parameter space points with 
50 < m x /GeV < 100 and 0.03 < p < 0.1. On the other hand, it is hard to satisfy our 
constraints in the region m x < 50 GeV, because the x 2 m this region rises sharply due to 
the Higgs boson mass and relic density constraints. The vertical branch in the figure is 
mainly due to the Higgs resonance effect, which can efficiently enhance the dark matter 
annihilation cross section when 2m x falls near m#. Hence, the coupling p has to be small 
correspondingly, in order to be consistent with WMAP data. On the other hand, when 



m x > mw, the XX ~^ W + W~ channel dominates the annihilation cross section [21. l3l. |29|. 
Therefore, we can see from the figure that in the 1 and 2 o C.L. bands of the horizontal 
branch the allowed p is roughly proportional to w? x (see Eq. (l34"j) at the Appendix). 

In Fig. [21 we show the profile likelihood on m x - 0^(0) panel against the experimental 
90% C.L. upper limit from XENON100. Clearly, the XENON100 data is only able to rule 
out 50 GeV < m x < 100 GeV. Current DM direct detection cannot constrain most of the 
parameters. On the other hand, the Higgs resonance region and most of the horizontal band 
can be tested in the future by XENON-IT (see the dashed line in Fig. [2]). 

Other than the Higgs resonance region, the WMAP constraint dominates the likelihood 
function as shown in Fig. ([2]), and therefore the largest likelihood of XENON100 only 
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FIG. 2. The profile likelihood of the spin-independent cross section <r™(0) for the SZ model 



projected onto the m x axis. The latest XENON100 limits |9|] are overlaid for comparison. The 



projected XENON-IT sensitivity is also shown. 



occurs at s <C b. Nevertheless, it is easier to satisfy the relic density constraint in the Higgs 
resonance region, and therefore the largest likelihood of XENON100 in the Higgs resonance 
region occurs at s — 1.0 such that s+b = o, by fine-tuning m x , p, and ran- As a consequence, 
the best fit of our scan appears in the Higgs resonance region. 



IV. PHENOMENOLOGY 

With the result of the likelihood determined, we can proceed to evaluate other observ- 
ables as predictions for the model, including gamma-ray lines, collider signatures, and muon 
anomalous magnetic dipole moment. 
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FIG. 3. The annihilation cross sections for the gamma-ray line from XX 77 (left) and XX lZ 
(right). 

A. Gamma-Ray Lines 



In Fig. 02]), we plot the cross sections for the gamma-ray line in the SZ model versus the 
profile likelihood projected onto the m x axis. The left panel is for XX ~^ 77 while the right 
one is for \X ~~ The Fermi-LAT data associated with different halo profiles are also 
shown for comparisons. It is clear to see that the prediction for the XX ~^ 77 annihilation 
cross section allowed by the profile likelihood is well below the Fermi-LAT data while that of 
XX ~~ >* 1% is even further below the Fermi-LAT data. Hopefully, future better measurements 
made by Fermi-LAT can put a dent in the allowed profile likelihood. 



B. Collider Signatures 

If the invisible mode of H — > XX opens up, we should study its impact on Higgs search 
at the LHC; in particular its effect on the branching ratios of H — > 77, H — > WW* and 



ZZ*, which apparent 



and ATLAS data 
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y show some excesses over the background. Since the current CMS 



131 1 showed that the excesses seen in 77, WW*, and ZZ* channels are 

n 

consistent with the expectation of the SM Higgs boson of 125 GeV, □ we cannot allow the 
5 The WW* and ZZ* decay modes are slightly below while the 77 mode is somewhat higher than the SM 
predictions. 
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invisible decay mode to be too large; otherwise the visible mode would become inconsistent 
with the current data. 

It is easy to show that the branching ratio for a visible mode would be its SM branching 
ratio multiplied by (1 — B- mv ) where B inv is the invisible branching ratio defined earlier as 
Tinv/ (rsM + Tinv)- in our scan in the previous section, we had required the invisible branching 
ratio -Bi nv < 0.4 such that each visible mode is reduced by an amount less than 40% so as 
not to upset the current data. If the dark matter mass m x > m#/2, the Higgs boson simply 
behaves like the SM Higgs boson. 

From our scan result in Fig. [TJ a few typical points which have the likelihood within 68% 
C.L. (1-cr band) can be identified as follows: 

1. Point A: m x = 53 GeV, p = 0.02, rrin = 125.3 GeV. The invisible branching ratio is 
right at 0.4. The significance of this point is that the Higgs boson still has a large 
branching ratio into XX- The collider signature that we will discuss below consists of 
a large missing energy. This is the point with the maximum likelihood, shown by the 
star in Fig. [TJ 

2. Point B: m x = 84.0 GeV, p = 0.042, uih = 125.2 GeV. This point gives x a mass close 
to my/ and hence above the Higgs decay threshold. This is the point at the low end 
of the second branch. 

3. Point C: m x = 608.3 GeV, p = 0.189, m# = 125.3 GeV. This point gives a heavy x 
that is still consistent with direct detection limits. 

The most common search modes so far for the dark matter are the monojet and monopho- 
ton plus missing energies. In this model, monojet or monophoton production must go 
through the Higgs boson H, so that the only sizable production cross sections have to go 
via 

bb^ H -> XX , 99^ XX, 

in which we can attach a gluon to the b or g leg, or attach a photon line to the b leg. Since the 
6-parton luminosity is small and gluon-fusion is a loop process, the monojet or monophoton 
rate would be relatively small. 

Since the DM candidate x only couples to the SM particles via the Higgs boson, the x 
will preferably couple to the heaviest fermion. At hadronic colliders, one of the interesting 
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FIG. 4. A Feynman diagram showing mono-6 jet production with missing energies, 
processes is 

gb bHW b X X , (25) 

where the superscript (*) on the Higgs boson denotes that the Higgs boson could be on- or 
off-shell depending on the mass of x- Obviously, it is dominated by on-shell Higgs boson for 
inn > 2m x . If mn > 2m x and mu is lighter than 2mw, the Higgs boson will dominantly 
decay into a pair of x- The corresponding collider signature would be a mono-6 jet plus 
missing energies. The other possible signatures would be the associated production of the 
Higgs boson with a gauge boson W or Z: 

pp -»■ W (Z)H^ -»■ Iv (£ + r )xx ■ (26) 

The final state in this case would consist of a charged lepton or a pair of charged leptons 
plus missing energies. 

We calculate the event rates of mono-6 jet, single charged lepton, and a pair of charged 
leptons plus missing energies at the LHC-7, LHC-8, and LHC-14. We impose the following 
selection cuts for the b jet or charged leptons and the transverse missing energy 

p Tb > 30 GeV , \r] b \ < 2 ; p Te > 25 GeV , \y t \ < 2 ; i> T > 50 GeV . (27) 

The cross sections for the mono-6 jet, single or a pair of charged lepton plus missing ener- 
gies are tabulated in Table HI The largest cross section comes from mono-6 jet production. 
However, when we apply the £>t > 50 GeV cut the cross section mono-6 goes down 50 
times. After further imposing the 5-tagging, the event rate would only be handful. An- 
other interesting signature is the single charged lepton plus missing energies. Counting both 
negatively- and positively-charged leptons the cross section could be as high as 16 fb at 
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the LHC-8. Given the LHC-8 can accumulate 20 fb _1 each experiment, it would be more 
than 300 events each experiment. The ZH production would give, on the other hand, two 
charged lepton plus missing energies with a few times smaller event rates. 



TABLE I. Cross sections for mono-6 jet, single charged lepton or a pair of charged leptons plus 
missing energies arise from Higgs boson production followed by H — > XX- We used the point A 
(rn x = 53 GeV, p = 0.02, and ran = 125.3 GeV, B mv (H — > xx) = OA). The selection cuts are 
defined in Eq. (127]) . 

Cross sections (fb) 



Subprocess LHC-7 LHC-8 LHC-14 

gb^bH ->• bxx 4.6 6.3 10.4 

ud -»• W + H -> l + u X X 9.2 10 19 

du^W'H ^ £~i?xx 4.7 5.8 12 

qq^ ZH -> 1 + 1~XX 2.2 2.6 4.9 



Note that for other typical points of the model, e.g., points B and C, the invisible decay 
mode of the Higgs boson is closed, and therefore the decay is similar to the SM Higgs boson. 
The process in Eq. fl25l) will then give rise to 36 or bWW* final states, depending on the 
Higgs boson mass. The processes in Eq. (f26|) will give one or two charged leptons plus either 
hi or WW*. 

The SM background for the mono-6 jet plus missing energy would be similar to the 
current monojet search in ATLAS 30j and CMS 3j|, but now with a 5-tag on the monojet. 



The largest background 



30, 



3jJ comes from Z + j — > vv + j and W + j — > iv + j with 



minor contributions from ti, single top production, and QCD multijets when leptons or 
extra jets get missing down the beam. On the other hand, background events with single or 
double charged leptons plus large missing energy comes from WZ — > Ivvv or ZZ — > livv 
with minor contributions from ti and single top production. Precise estimations of these 
backgrounds are beyond the scope of the present paper. 



1(3 




FIG. 5. Two loop Feynman diagram contributing to the muon anomalous magnetic dipole moment. 
C. Muon Anomalous Magnetic Dipole Moment 

The experimental value of the muon anomalous moment a M = (g M — 2)/2 is 

a° xp = 116 592 089(63) x 10~ n , (28) 

while the SM prediction is 

a™ = 116 591 802(49) x 10~ n . (29) 

The 3.6a discrepancy between the above experimental measurement and theoretical calcu- 
lations based on using the e + e~ annihilation cross section for the estimation of the hadronic 
correction 32] 

Acv = a^ xp - af 1 = 287(80) x 10~ n (30) 

could be a harbinger of various new physics beyond the SM. The contribution to the muon 
anomalous magnetic dipole moment in the SZ model first shows up at the two loop level 
(See Fig. [5]). Detailed expressions can be found in the Appendix. In Table (JTTJ) , we show 
the numerical results of a M for the three typical points A, B and C from our scan. For all 
the relevant parameter space, we have checked that the contribution is negative and many 
orders of magnitude below the current experimental sensitivity. 



V. CONCLUSIONS 



The simplest dark matter model is realized by 
model as was discussed quite a while time ago in 



iaaing a real scalar singlet to the standard 
l|, long before the popular dark matter 
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TABLE II. Muon anomalous magnetic dipole moment for Points A, B and C of the likelihood. 



Muon Anomalous Magnetic Dipole Moment (a^) 
Point A Point B Point C 

-1.47 x 1(T 21 -3.19 x 1(T 21 -1.67 x 1(T 21 



candidate of neutralino in MSSM model took the central stage. In this work, we use the most 
current experimental constraints of the relic density from the 7 year WMAP data, latest 
XENON100 data, annihilation cross sections from Fermi-LAT based on 10 dwarf spheroidal 
satellite galaxies of the Milky Way, as well as the 125 GeV Higgs candidate as discovered 
recently by the LHC, to pin down the profile likelihood for the parameters p and m x of the 
model. 

The collected points are then used to evaluate the cross sections for the gamma-ray lines 
from XX ~ >* 77 an d lZ and found that they are well under the current limits from Fermi-LAT 
data. A small part of the allowed parameter space around m x ps 70 GeV barely touches the 



Fermi-LAT data with Einasto halo profile. Recently, an interesting analysis in Ref. : 



using 



the Fermi-LAT data suggests there could be a gamma-ray line around 130 GeV that may be 



related to dark matter annihilation. However, other authors 
sources like the fermi-bubbles 



34j suggest that astrophysical 



351 ] could also be responsible for this line signal. The gamma- 
ray lines in this simplest dark matter model cannot accommodate this line signal based 
on the profile likelihood determined by the global fitting with the experimental constraints 
mentioned above. 

We also study the LHC signals of mono-6 jet, single charged lepton or a pair of charged 
leptons plus missing energies of the model. The most interesting case is the single charged 
lepton plus missing energies which can arise from associated production of WH followed by 
W — > Iv and invisible decay of the Higgs. With a luminosity of 20 fb _1 for each experiment 
of ATLAS and CMS at LHC-8, we expect several hundreds of such events based on the Point 
A. 

We also evaluate the muon anomalous magnetic dipole moment of the model and found 
that it is many orders of magnitude below the current experimental limit for all relevant 
parameter space. 

More stringent constraints are expected for this simple model of dark matter as more 
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data from the LHC, direct and indirect detection experiments become available in the near 
future. 



APPENDIX 



1. Matrix Elements 



In this Appendix, we list the matrix elements and annihilation cross sections for all the 
two body processes needed in the calculations of the relic density and indirect detection. 
Let s to be the center of mass energy given by s = 4m^/(l — where v x = 2f3 x with /3 X 
being the velocity of the dark matter. No is the color factor, 1 for leptons and 3 for quarks. 
(1) XX -+ ff: 



2N c p 2 m 2 f 



spin 



(2) XX -»• VV (V = W or Z): 



m 



Amj 



2 \ 2 , 2 t->2 \ 1 

m H ) +m H T H j 



(31) 
(32) 



spin 



(s 



m 



1 



(TV, 



1 + SyZ 167T 



Amy 



2 \2 
H) 

1 

2 



1-4 



mi 



12 



mi 



m 



2 \2 
HI 



m\Y\ 



m 



v 



+ 12 



mi 



(33) 



(34) 



Here 5yz is a Kronecker delta to account for the Bose statistics of the ZZ final state. We note 
that MicrOMEGAs computes process cross sections by CalcHEP 36(. However, we found 
that the amplitude squared of XX ~^ W + W~ /ZZ differs between Eq. ( 1331) and the result 
from CalcHEP. The factor inside the square bracket of Eq. (|33|) is 



-4 



1-4| 
f 12 



while the corresponding factor in CalcHEP reads 
discrepancy we rescale the cross section by the ratio of these two factors 
(3) XX -> HH 



Due to this 



spin 



3m 2 H 



pv 



pv 4 



[s — m 2 H ) + im H T H t 



mi 



u 



mi 



(35) 
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(4) XX ->• lZ 



647TS 



1 - 



Am 2 H 



dcos0^|A<| : 



spin 



p 2 v 2 



is — m 



2\ 2 



■\Ay Z (s) 



spin 
2 



lQn 2 m w 



2 (s - m 2 H f + m 2 H T 2 H 

^-2Q f (T'f L - 2Q f sin 2 9 W ) f 

4 cos e w r w + V VM 7 — ^ ^iv£/ 



(36) 



(37) 
(38) 



i (3_ tan 2 ^)m 2 1/ + x?/r(-5 + tan 2 ^)m^-i(l-tan 2 ^)sl 
ax / ay 



m 2 w - y(l - y)m 2 z + xy{m 2 z - s) - i0+ 
(3 - tan 2 9 W )I (r w , r zw ) + [(-5 + tan 2 6 W ) - 2(1 - tan 2 6 w )t w ] J(t w , t zw ) (39) 



1 pl—X 

dx 

Jo 



dy- 



(4xy — 1) m 2 . 



m 



/ -2/(1- y)m% + xy( 

4J(Tf,T Z f) - I(Tf,T Z f) 



m z -s) 



i0+ 



(40) 



Here, ry^ = s^m 2 ^, r/ = s/4m 2 , r^w/ = m z /^ m wi an d = ^|/4m 2 . I and J are given 
by 

I(rnT 2 )= 1 -(T 1 -T 2 y 1 (f(T 1 )-f(r 2 )) (41) 



•tyi, r 2 ) = -1 (n - r,)- 1 [1 - (n - r 2 )- x (/ ( Tl ) - / (r 2 ))] 



1 _ 2 

+ ^t 2 (ti-t 2 ) (»(ri)-^(r 2 )) 



with / and g denned by 



f(r) 



[sin 1 y/r\ ' 



In 



( i+VT- 



\i-VT- 



for r < 1 



for r > 1 



g{ 



T 



1=^ tan" 1 



l-r 



<71\ 



1 / m|\ 5 p 2 t) 2 (s — m|) 2 



for r < 1 
for r > 1 



|Ayz(s) 



(5) XX -»■ 77 



(42) 



(43) 



(44) 



(45) 
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E l-^l 2 = P ^f- ( 2 1 2r2 IA^)! 2 (46) 



spin 



V = 2 + 3r^ + 37V 1 (2 - r^ 1 ) / (r w ) (48) 

7; = -2TJ 1 [1 + (1 - Tj 1 ) f (T f )] (49) 

with / defined in Eq. (l43|) . 



°v x = -±-pv- — ; 2r2 \^)\ 2 (so) 

o47T (s — rriff) + m z H T z H 



(6) XX -> 55 



2^ \ M \ = 7 ^ ^I^M 5 )! ( 51 ) 



spin 



-P 2 « 2 7 77^ ^K».»M| 2 (53) 



9 

where J g is given by Eq. (|49l . 

1 (N 2 C -1) 



X 64vr 4 (s - mf,) 3 + m^rf, 



2. Muon Anomalous Magnetic Dipole Moment 



Following a similar procedure as in the QED case 37J, one can readily obtain the following 
result for the muon anomalous magnetic dipole moment 

^ = f «« s / ^ f d2 l0E [1 + x[z)e] (54) 

where 

H(x)=x 2 + (l-x)r 2 Hfl , (55) 

*(z) = z{l - z)rl x , (56) 

with r 2 H = Tnjj/rn 2 and r^ x = m 2 /m 2 ^. Performing the £ integral, we end up with a 
two-dimensional integration for 

a u = — - [ dx [ dz n , ' . (57) 

128Wo 7 -D(x, z) 1 ; 
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where 



N(x, z) = x 2 (l -xf 



and 



2(l-H(x) X (z)) 
H\x) X \z) 



D(x, z) = H 2 (x) 



1 



H(x)x(z) 



H(x)x(z) 



- 1 



log (H(x)x(z)) 



(5* 



(59) 
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